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The lamellar phase in diblock copolymer systems appears as a result of a competition between 
molecular and entropic forces which select a preferred periodicity of the lamellae. Grain boundaries 
are formed when two grains of different orientations meet. We investigate the case where the lamellae 
meet symmetrically with respect to the interface. The form of the interface strongly depends on 
the angle, 6, between the normals of the grains. When this angle is small, the lamellae transform 
smoothly from one orientation to the other, creating the chevron morphology. As 9 increases, 
a gradual transition is observed to an omega morphology characterized by a protrusion of the 
lamellae along the interface between the two phases. We present a theoretical approach to find 
these tilt boundaries in two-dimensional systems, based on a Ginzburg-Landau expansion of the 
free energy which describes the appearance of lamellae. Close to the tips at which lamellae from 
different grains meet, these lamellae are distorted. To find this distortion for small angles, we use 
a phase variation ansatz in which one assumes that the wave vector of the bulk lamellar phase 
depends on the distance from the interface. Minimization of the free energy gives an expression 
for the order parameter 4>{x,y). The results describe the chevron morphology very well. For larger 
angles, a different approach is used. We linearize <f> around its bulk value <f)L and expand the free 
energy to second order in their difference. Minimization of the free energy results in a linear fourth 
order differential equation for the distortion field, with proper constraints, similar to the Mathieu 
equation. The calculated monomer profile and line tension agree qualitatively with transmission 
electron microscope experiments, and with full numerical solution of the same problem. 
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I. INTRODUCTION 



The lamellar phase is one of the possible phases with spatial modulations that can be found in a wide variety 
of physical and chemical systems. These include diblock copolymer melts, mixtures of diblock and homopolymers, 
aqueous solutions of lipids or surfactants, Langmuir monolayers, and magnetic garnet films 0. Modulated phases are 
the result of a competition between forces, one of which prefers ordering characterized by a non-zero wavenumber, 
while the other prefers a homogeneous (disordered) state. Below we shall employ the language appropriate to block- 
copolymers, but our work applies equally to other systems. 

We consider diblock copolymer melts in which the two polymer blocks are incompatible. This incompatibility is 
characterized by a positive Flory parameter \- Because of the covalent chemical bond between the A and B blocks, 
the system cannot undergo a true macrophase separation. Instead it undergoes a microphase separation characterized 
by A- and B-rich domains of a finite size. Various modulated phases such as lamellar, hexagonal, and cubic are 
observed [§-|tJ depending upon the % parameter as well as on the relative lengths of A and B blocks. In most cases, 
the x parameter depends inversely on the temperature. Hence, the lamellar and other modulated phases will not be 
stable at high temperatures and the polymer melt will be in a disordered state. 

The equilibrium behavior of diblock copolymers in the bulk is by now well understood |3]-f7) ■ Single domain bulk 
phases, however, are rarely observed in experiments because it is extremely difficult to completely anneal defects. 
In most cases, due to very slow dynamics and energy barriers, the microstructural ordering is limited to finite-size 
domains (or grains) separated by grain boundaries. These defects are very common to block copolymer systems and 
are readily observed in experiments HQ ■ 

Because these domain boundaries and defects are so abundant in polymer melts, it is of interest to study their 
energetics and other characteristics. In this paper we concentrate on the relatively simple situation of domain bound- 
aries in lamellar phases in which there is no twist between the two grains, only a tilt, as is shown schematically in 
Fig. la. The system is translationally invariant along the z-direction and can be described by its x-y cross-section 
only, reducing it to an effective two dimensional system. As can be seen in Fig. la, the distance between lamellae 
along the grain boundary, d x , is larger than the lamellae spacing in bulk d, by a factor of l/cos(0/2). This causes an 
increase in the local free energy density related to the grain boundary. 

The experiments of Gido and Thomas H and Hashimoto and coworkers B show that the response of the system to 
this increase in free energy depends strongly on the tilt angle 9 between the grains. The so-called chevron morphology 
H occurs when the angle is small. In such a situation, lamellae transform smoothly from one orientation to the 
other, creating V-shaped tips. The rounding of these tips reduces the interfacial area between lamellae at the expense 
of introducing a curvature energy. For block co-polymers the profile shape at the tip is determined by the local 
relaxation of the stretched chains. The chevron morphology is shown (schematically) in Fig. lb and again in Fig. 2a. 
A gradual transition to an omega-shaped tip is observed when 9 is increased, as seen in Fig. 2b. This morphology 
is characterized by protrusions of the lamellae along the interface between the two phases. The protrusions can be 
understood as a different attempt of the system to reduce the cost of the boundary, while still complying with the 
geometrical constraints. Essentially the system tries to create a lamella similar to those in the bulk but which is 
aligned along the interface itself. 

The basic phenomenology of these grain boundaries was presented by Gido and Thomas |J. Netz, Andelman, and 
Schick || then considered the phenomenon employing a Ginzburg-Landau free energy which was minimized numer- 
ically, and obtained both the chevron and omega morphologies. Matsen |8[ considered the block-copolymer system 
explicitly and employed self-consistent field theory. Not only did this produce the chevron and omega morphologies 
for small and intermediate angles, respectively, but also a symmetry-broken omega for large enough angles. In the 
latter case, the response (and shape) of lamellae composed of one of the blocks differs from that of the other. Such 
symmetry broken boundaries were indeed observed by Gido and Thomas H . 

In Sec. H we adopt the Ginzburg-Landau free energy functional employed earlier [Q]. The advantages to using this 
functional are its simplicity and generality, while retaining the essential ingredients that capture the behavior of the 
system. In contrast to the complete minimization of the free energy functional which requires a numerical calculation, 
we shall employ here a simple ansatz for the form of the grain boundary in order to obtain analytic results. Recently, 
similar methods were employed to obtain analytically the interface between the lamellar and disordered phases of 
diblock copolymers || . Our motivation is to demonstrate that the essence of these interesting morphologies does not 
depend on strong segregation conditions or a large number of Fourier components, and so should be observable in all 
systems with modulated phases. 

The chevron structure is obtained by using the bulk lamellar phase solution with a constant amplitude, but with 
varying wavevector. Minimizing the free energy subject to the proper geometrical constraints, we find an equation for 
the wavevector. This is done in Sec. III. Beyond the chevron regime (i.e. for large inter-domain angles), this approach 
will not be adequate, because the amplitude of modulations will have to vary as well. To this end we expand <fi(x, y) 
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around the two bulk lamellar phases, giving rise to an equation for a small distortion field. Close to the interface, 
the sharp tips of the lamellar phase (see Fig. la) are smoothed-out and the protrusion characteristic of the omega 
morphology appears. Far away from the interface, the disturbance vanishes and the bulk lamellar phase is recovered, 
as is shown in Sec. 



IV 



In Sec. IV A we discuss some features of our method, the analogy to, and the differences from 
the Schrodingcr equation for electrons in a one-dimensional periodic potential (known also as the Mathieu equation). 
In Sec. [v] we report our results, and discuss them in Sec. [vj. 



II. THE MODEL 



An order parameter ip(R) — [i/)a(R) — V'.b(R)] 1S defined as the difference in local A and B monomer volume 
fractions. We employ the following Ginzburg-Landau free energy functional of this order parameter: 
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where ks is the Boltzmann constant, and T is the temperature. With c\ and ci positive, the first two terms favor 
a uniform, disordered state. The coefficient of the third term, C3, is negative and therefore induces the system to a 
modulated, ordered phase. The Laplacian squared term ensures that these modulations are not too large. This type 
of free energy functional (and some variants of it) has been successfully used to describe the bulk phase of diblock 
copolymers |5]||, amphiphilic systems |1C| ], Langmuir films and magnetic (garnet) films Jt2[ . 

Assuming an infinite system that is non-uniform only along one direction, one can minimize this free energy to 
obtain a solution describing the lamellar phase ipL oc cos(Q • R). One readily finds that the transition to the lamellar 
phase first occurs at a wavenumber Q = (— C3/2C4) 1 / 2 . It is convenient, therefore, to introduce the dimensionless 
position vector r via 
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so that the dimensionless and rescaled free energy functional takes the form 
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The coefficient ci changes sign as the interaction strength increases or temperature decreases. In a polymer system, 
ci can be shown to be proportional to 1 — Xf/x*fi where xj is the Flory parameter |l3| ], and x*f the value of this 
parameter at which the coefficient c\ passes through zero [ ft4f . 

To make contact with Ref. O, we set the ratio of coefficients C1C4/C3 to be equal to 1 — x- For this particular choice 
of coefficients, the free energy functional becomes 



with 
2^ 



The bulk lamellar phase, with wavevector in the x direction, is described by 

4>l = 0o + 4> q cos(qx) 
the average volume fraction difference 



(6) 
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q — 3/4 the amplitude of the variations 
namical stable phase in a range of the system parameters: 
difference in the average volume fractions of A and B monomers 



= Q(— C3/C4) 1 / 2 = 1/V2 the optimal wavevector, and <f> q = 
It can be shown l^-Q that this lamellar phase is the thermody- 
'0 and x- F° r block copolymers, 4>o is proportional to the 
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Upon substitution of 0^ into Eq. ^[ one obtains the free energy per unit volume of the bulk lamellar phase 
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where V = J d 3 r is the rescaled, dimensionless, volume of the system. 

Some remarks are now in order. First, as was mentioned in the introduction, the free energy (|^) can have other 
non-lamellar modulated solutions P,|l5||. We will not consider them in this paper since our aim is to study defects 
inside lamellar phases. Second, the validity of a single optimal mode can be justified in the weak segregation limit 
(i.e., near a critical point or a weak first-order transition). Far from the critical point higher harmonics are needed to 
describe the optimal lamellar phase [§,1). In addition, very close to the critical point corrections due to fluctuations 
are important |16 i*7f|. 



We now turn to the tilt-boundary problem, where two lamellar domains, both lying parallel to the x-y plane, meet 
with an angle 9 between their normals. (See Fig. la.) The x-axis is along the line interface between the two lamellar 
domains. The y-axis is perpendicular to it. In these variables the lamellae in the two grains are described by 



0L = 00 + 9 cos(q x x =F q y y) 



(9) 



where q x = qcos (0/2) and q y = qsm (9/2) are the components of the optimal wavevector, q = (q x , q y ). Their inverses 
provide characteristic length scales in the x and y directions, respectively. The y < half plane is a reflection through 
the x axis of the y > half plane, so it is sufficient to consider only the upper half plane, y > 0. The system is 
periodic along the x axis, with wavelength d x = 2ir/q x = 2ir/[qcos(9/2)]. 



III. DESCRIPTION OF THE CHEVRONS 



For small tilt angles, the lamellae transform smoothly from one orientation to the other, showing the chevron 
morphology. We will assume that for small enough tilt angles, the only change of the functional form of the order 
parameter is through the wavevector p8|,^9[ . The aim of this section is to show that, in the chevron regime, the 
behavior of the system close to the interface is, in essence, quite similar to that far from the interface. The diagonal 
lines in Fig. la show each bulk lamellar phase in its respective half-plane, and the sharp tips which result from their 
intersection. These sharp tips will be smoothed out in the chevron morphology, Fig. lb. 

We use the following ansatz for the order parameter: 

c (r) = 0o + g cos[q x x + q x u(y)}, (10) 

where the amplitude 4> q is identical to that in the bulk solution 0^. This choice is motivated by the fact that, in the 
chevron morphology, the amplitude of the order parameter appear to be rather constant, while it is the phase which 
changes smoothly from one grain to the other. The local direction and magnitude of the wavevector depends on the 
distance y from the s-axis. Far away from the interface, the lamellae must return to their bulk orientation, implying 
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and at the interface, the continuous function u(y) satisfies 

lim u(y) = 

Symmetry with respect to inversion across the a;-axis means that 
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= -<f> q q x u'(0)siii(q x x) = 



which implies 
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We insert the form ( |io| ) into the free energy functional (^). The integration over x and z can be carried out. 
The remaining integration over y shows that F[0 C ] is proportional to the dimensionless, rescaled, volume V of the 
system. This reflects the fact that the order parameter profile approaches its bulk value far from the grain boundary. 
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Subtraction of the bulk free energy produces a functional for the free energy per unit area of the grain boundary. 
After simple manipulation, one finds the expression 
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From the boundary condition, Eq. |Tl], one sees that far from the grain boundary s attains the value ±1 and that its 
derivative vanishes. It is clear, therefore, that if s approaches its limiting value sufficiently quickly, the grain boundary 
free energy is finite. This is indeed the case. 

The Euler-Lagrange equation which minimizes the grain boundary free energy is 



^— + s (l- s 2 )sin 2 (0/2) = 
dy 2 
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which has the solution 
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and 



u (y) = -f — tan (9/2) In cosh (q y y) 
1y 

We have thus found the order parameter profile following the initial ansatz ( Jlp) ) 

<j> c {x, y) =<f> + (j) q cos[q x x - In (2 cosh^y))] 
The grain boundary free energy can now be obtained by inserting the profile solution, Eq. (]2Cj) , into Eq. (fL5l 
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sin 3 (6>/2) 



(19) 
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The order parameter we have calculated describes the chevrons very well. The tips of the V-shape structure are 
rounded-off, and far away from the interface the bulk phase is restored. The values of the grain-boundary energy are 
close to those obtained from full numerical minimization of the free energy functional (pf) , see Ref. @. The expression 
( pl| ) for the grain boundary energy shows a 9 3 scaling for small angles Jl8| 

The width of the grain boundary is the characteristic distance over which the phase of the order parameter profile 
deviates from its bulk value. From the profile of Eq. |2(], we see that this distance is l/q y which, for small grain 
boundary angles, varies as 1/9, in accord with well established results p8) . 

The deficiency of the above approach is that it does not give the cross-over from the chevron to the omega mor- 
phology. For this end another approach will be used in the next section. 



IV. DESCRIBING THE OMEGAS 



As the tilt angle 9 is increased, the chevron structure is deformed more and more. For large angles the lamellae 
protrude along the interface between the grains, creating the Jl-shaped structure. In this regime the polymer profile 
is qualitatively altered, and calls for a different approach. 

In this section, we use the bulk-phase solution as a zeroth order approximation and determine a correction to it; 
that is we write 
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<p(x, y) = <p L {x, y) + 5cj)(x, y) 



(22) 



Of course our ansatz for the chevrons, which was motivated by the smooth variation of phase which they display, 
can also be written in this form with a particular choice of Sip. Below we shall choose a different Sip based on the 
observation that the distortion of the omegas is well localized near the grain boundary. 

After substitution of this form into Eq. (0) the free energy can be written as a sum of two parts: 



F[(P] =F L + 1 A 
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where Fl — F[<pz] is the bulk free energy, proportional to the volume and is the grain boundary energy proportional 
to the area of the boundary. To second order in Sip, the latter is 
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Since <pL minimizes Fl, we need to find the function 8<p{x,y) that minimizes 7. This, in principle, is done via the 
Euler-Lagrange equation. 

The boundary conditions for <p(x,y) follows from the symmetry of the grain boundaries, and from the requirement 
that ip approach its bulk value away from the interface: 
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(25) 
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In the above, n is odd. These conditions impose boundary conditions on Sep, because 4>l is known, and <p — <Pl + Sip. 

The distortion field S<p(x, y) will be found based on an ansatz. Let us evaluate the y-derivative of the bulk tilted 
lamellar phase: 
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Because the distortion of the omegas is well localized near the grain boundary at y = 0, it is therefore reasonable to 
assume that <p has the following form: 



<pn{x, y) = tp L (x, y) + f(y) sin (q x x) 



(28) 



The procedure we adopt is as follows. This ansatz for <p is inserted into the free energy density and the integration 
over x and z is carried out. The grain-boundary energy functional 7 then depends on the unknown amplitude /(y). 
This is minimized by an Euler-Lagrange equation, which results in a fourth order ordinary differential equation for 
f(y). The boundary conditions, Eqs. (p^-(Pq), will translate into boundary conditions on f(y). The solution of the 
Euler-Lagrange equation gives, in principle, everything we want to know about the system: spatial distribution of the 
order parameter, line tension, etc. 



We begin by putting our ansatz (|28|) into the expression 1 
yields 
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The Euler-Lagrange equation for the function f(y) is obtained by minimizing (E 
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[A + C cos (2q y y] ) / + Bf" + f"" = (30) 
where A, B and C depend on the parameters of the problem as follows: 

A{e) = i-x + <t>t + y q + ii-<& 

B(9) = 2ql 

Note that the profile equation (BOf) is linear in / because the free energy, Eq. fl24) is second order in S(f>. Equation 
( |30| ) is similar to the Mathieu equation, which is the Schrodinger equation for an electron in a periodic (sinusoidal) 
one-dimensional potential as appears in many solid state physics problems. The parameter A plays the role of the 
electron total energy, and C is the amplitude of the periodic potential. Unlike the solid-state case, here the "energy" 
parameter A, the "kinetic term" parameter B as well as the periodicity of the "potential" term depend on the angle 
9, leading to a more complex band structure as function of 9. Lastly, the differential equation is of fourth, not second 



order. Another useful observation is that equation (30) is invariant with respect to the spatial variable transformation 
y — ► —y. Therefore, the solutions can be classified as symmetric and asymmetric. As we will see, the symmetric and 
asymmetric solutions do not satisfy separately the boundary conditions, so a combination of them will be needed. 

Since this equation is linear with periodic coefficients, a solution to it will have the Bloch form. This is also known 
as the Floquet theorem pfj| : 

f(y)=e k y g (y) (31) 

where g (y) is a periodic function with period d v = it/q v which is the same periodicity as appears in Eq. (|30|). Hence 
it is appropriate to write it as a Fourier series: 

oo 

g(y)= £ a n e 2in **y (32) 
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Substituting the Bloch form Eq. ( |3l| ) into Eq. (|30"|), we get a sum of exponential terms. Demanding that the coefficients 
in front of every exponent vanish, we obtain the following recursion relation: 

(A + B (k + 2inq y f + (k + 2inq y f^ a n + (a„_i + a n+1 ) = (33) 

The appearance of a„_i and a„ +1 is due to the cos(2q y y) term in (|30|). At first glance, it seems that for every k, 
choosing "initial values" for the coefficients {a n } gives a valid solution. However, a closer inspection shows that for an 
arbitrary k vector the series a n will diverge. Only a very specific value of k (eigenvalue) will give a convergent series. 

The method by which we find this value is as follows (e.g., see Ref. pG]). Rewrite the recursion relations (B3h as 
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On-i A + B(k + 2inq y ) + (k + 2inq y ) +\Ca n+ \/a n 



for n > 0, and similarly 



in+i A + B (k + 2inq y ) 2 + (k + 2inq y ) i + hCa n -i/a n 



(34) 



(35) 



for n < 0. For very large values of n, a n +i/a n should be much smaller than one, so one can start from some N 3> 1, 
assuming that a^+i/aN 0, and get 

a ziZ ( 36 ) 

a N -i A + B{k + 2iNq y f + {k + 2iN y ) 4 

then iterating backward gives the ratios {a„/a„_i} for < n < N. Carrying out the same procedure for negative n's 
one arrives at the stopping condition: 

A + Bk 2 + k^--c(^ + ^-) (37) 
2 \ ao ao J 

from which k is deduced, because the ratios ai/ao and a_i/ao are already known from Eqs. (|3^)-(p5|). The iteration 
scheme is to choose an initial guess for k, put it in (E^), and use (EmI) and (|35|) successively for all n ^ 0. When n = 
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is reached, k is calculated from (J37|), and put back in (pq). The process repeats until convergence is achieved. From 
the required boundary conditions at infinity, only fc's such that Re(fc) < are acceptable, recalling that the Bloch 
form (|3l|) contains an e ky term. 

It can be seen from ([m]), ( ^5|) that if a certain value of k gives a physical converging solution for y — > oo, then fc*, 
the complex conjugate of fc, will also be a convergent solution, but with amplitudes a*_ n . Solutions with —k and — fc* 
are possible, too, but are discarded because for them f(y) diverges as y — > oo. The functions with definite symmetry 
include both k and ~k and hence diverge at infinity. Consequently, by dropping the — k solution we choose a specific 
combination of the symmetric and asymmetric functions with respect to y. The function / is a combination of the 
two independent solutions, and since it is real, it must be equal to: 



f(y) 



Mi 



E 
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+ C.C. 
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A choice of k such that Re(fc) < ensures that the disturbance will decay away far enough from the interface, 
y — > oo. The use of a linearization scheme leaves us with a linear ordinary differential equation, so that we lose the 
ability to impose all bound ary conditions, and forces us to use only the first and third derivatives of /, which by the 
use of Eqs. (@), © and (f§) are 



dy 
dy 3 



-q%<t> q 



(39) 
(40) 



In summary, the series {a n }) can be determined from the recursion relations (35), (|3J), while the only remaining 
unknown is the complex parameter ag. The function f(y) itself depends on this series as well as on the two boundary 
conditions. Substitution of ([38]) in the above boundary conditions fully determines both the phase and magnitude of 
ao and therefore sets the complete solution to the problem. 



A. Band-gaps and degeneracy 



In the usual Schrodinger equation for periodic potentials, one encounter energetic gaps in the energy spectrum. 
They occur whenever the fc-vector crosses the edges of a Brillouin zone. It should be expected that similar phenomenon 
happens here, too. Indeed, in our case it will happen whenever 

Im(fc) = mq y (41) 

for some integer m. For k fulfilling ( |4"l| ) there is a degeneracy in the two previously found solutions. It can be seen by 
noticing that if ( |4l"| ) holds, k and k* differ by 2imq y . Therefore, the amplitudes {a n } corresponding to the eigenvalue 
k* are the same as those {a„_i} corresponding to k, and the two solutions are dependent. If m is even, the two 
solutions are real. One can see this by using k = k r + 2inq y (n being an integer) and noticing that 



c (k r +2inq y )y 



J2 a n e 2mq y y = e krV a n _ l e 2inq ^ v (42) 



Using the right term of ( |42| ) in the iteration process it is clear that a n = (a_„)*, and f(y) is real. 

In these band gaps there is a splitting of the real part of k; another independent fc-vector appears whose imaginary 
part is the same, but which has a different real part. This is clearly seen in Figure 8a, where the band gaps appear 
to be centered around 18°, 26° and 42°. The last gap occurs for 9 > 127°. 

The two solutions constructed by the Bloch form |3l]) are obviously independent. It should be noted that the 
appearance of energy gaps is a mathematical consequence which does not introduce any singularity into the physical 
grain-boundary energy. It is due to our approximative linearization scheme, and is not expected to occur if non- 
linearity would have been included. 
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V. CHEVRONS AND OMEGA-LIKE PROFILES 



We discuss the results obtained by using the analytical eigenvalue scheme detailed above. For small angles (small 
tilt) the response of the system is weak. Namely, the lamellae gradually change their orientation and the profile is 
well described by the analytical chevron form, Eq. (^0|). At higher tilt angles the polymer chains at the interface 
are stretched much more than their preferred length. Here our analytical scheme (Sec. IV) gives rise to profile shape 
quite different than those obtained in Sec. Ill for the chevron. 

Figure 2 shows contour plots of the order parameter profile (f>n(x,y). We have taken the average of the order 
parameter, <fi , to be zero. The gray levels indicate the magnitude of the order parameter: black regions corresponding 
to the maximum values of <pQ (rich in A polymer), while white to its minimum value (rich in B polymer). The 
interaction is set to \ — 1- I n Fig. 2a the chevron morphology is represented using the results of Sec. IV, with tilt 
angle 8 — 20°. The plot is identical when using the expressions of Sec. III. A smooth changeover between the two 
lamellar phases is observed. The chevron solution follows nicely the analytical form of Eq. (po|). For large angles 
(8 = 130°), the omega structure takes over as is shown in Fig. 2b, with large protrusions of the lamellae at the 
interface. As is explained below, for large angles the excessive packing frustration of the chevrons cost more energy 
than the omegas. Chevron and omega morphologies are also shown in Fig. 3 but with x = 0-76 much closer to its 
critical value of 3/4. As one goes away from the y = interface, undulations of <pQ are encountered. There are more 
of them as the order-to-disorder transition (ODT) is approached, and the omega structure is more evident. 

Figures 4 and 5 show the same plots, but with the lines of interface enhanced. The regions in white, gray and black 
correspond to <pn < —0.2, —0.2 < (f>n < 0.2, and <f>a > 0.2, respectively. The gray marks the interface between the A 
and B rich regions. The equi-^o lines clearly show the form of the interface in the chevron and omega morphologies. 
Notice that as the ODT is approached (\ — 0.76 in Fig. 5), there are more undulations apparent on top of the bulk 
lamellar phase. In the chevron regime, the A/B interfacial width remains almost uniform, while in the omega regime 
it varies close to the kink ||. Fig. 4 shows quite clearly that the leading and trailing edges of the lamellae are quite 
different when the angle of the grain boundary is large. This is emphasized in Fig. 6, for 8 — 130° and x — 1-0, the 
same conditions of the lower plot of Fig. 4. The lines indicate where <fi = 0.5. That the protrusion at the leading edge 
is much more pronounced than at the trailing edge is quite reminiscent of the profiles seen in experiment Q with one 
exception; in experiment, only half of the lamellae look this way, the other half hardly display protrusions at all; that 
is, the symmetry between positive and negative order parameter domains is broken. 

Results for the grain-boundary energy 7tb are shown in Fig. 7. The x parameter is arbitrarily fixed to be x — 0.76. 
We show the grain-boundary energy calculated analytically by means of the methods of Sec 



III 



and which is valid 



for small 8, as well as that calculated in Sec. |fV| The grain-boundary energy ttb is an increasing function of the 
angle 8. In the small angle regime the grain-boundary energy obtained in Sec. [II scales as 7 ~ 8 3 . This scaling is 



also satisfied (to a good approximation) by the solutions obtained in Sec. IV. To see this, the same data is plotted in 
a different fashion in Fig. 8. In 8a, 7 1 / 3 is plotted as function of 9, while in 8b the data is plotted on a log-log plot. 



From both parts of Fig. 8 we conclude that the solution obtained in Sec. IV gives a power law with exponent of 2.91. 
Accuracy of the leftmost point of the solid curve in Fig. 8b is doubtful, due to poor convergence of the numerical 
iteration scheme for small angles, 8 < 5°. As the tilt angle grows, deviations from the 8 3 behavior become larger, and 
the omega morphology, with its lower energy, appears gradually. 

For intermediate and large angles, our results are supported by a full numerical solution of this problem [Q. 
Presumably, for small angles there is agreement too. This needs to be further checked since in Ref. || the smallest 
angle was 8 = 28°. Note that in contrast to the full numerical solution in which the order parameter profile was 
obtained via a functional minimization, here we employed numerical means only to obtain the value of the eigenvalue 
k, while the profile equation was solved analytically. 

In Fig. 9 we show the A;-vectors found by the use of the iteration scheme, as a function of tilt angle 8, with x — 1 
fixed. Outside of the band-gaps, k and k* are valid solutions, so for clarity only the k with Im(fc) > is shown. Notice 
that Re(fc) is always negative. The band-gaps are clearly seen on this graph as regions where Re(fc) have two distinct 
values (hollow circles). A check on the imaginary part of k in these regions reveals that it is an integer multiple of q y . 

Our results agree well with experiment with the exception that we do not obtain the symmetry-breaking transition 
of the omega morphology ||. To do so, one must add to an ansatz for the order parameter at least a term which 
varies as 2q x x, in addition to the fundamental term varying as q x x. 

VI. CONCLUDING REMARKS 

We have used a simple Ginzburg-Landau free energy functional to investigate the profiles between lamellar phases 
of diblock copolymer. Our analytic results give good qualitative agreement with experiment in the weak segregation 
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regime, and with full numerical solution of the same free energy model. The observed chevron morphology develops 
gradually into an omega morphology for intermediate tilt angles. For small angles, the use of a periodic order 
parameter with constant amplitude but varying wavevector suffices to describe the order parameter profile. This is 
well described by the chevron morphology. For intermediate angles the change of the profile at the interface deviates 
significantly from the bulk, and requires a different treatment. The deviation from the bulk lamellae was found, 
and gave rise to the protrusion characterizing the omega morphology. The symmetry breaking of this phase was not 
obtained. 

We were able to calculate grain-boundary energies and to determine that they scale as the cube of the angle [p||l§|l 
for small angles. As the tilt angle grows into the intermediate regime, the profiles deviates continuously from the 
chevron shape. As the angle approaches 180°, the energy must go to zero because the grain boundary vanishes at 
that angle. This does not occur in our analytic derivation, since, for such large angles, the linearization assumption 
is no longer valid. 

Interesting extensions of the present work would describe interfaces of two perpendicular lamellar phases (so-called 
T-junctions) , interfaces between modulated phases of other symmetries, and the inclusion of a twist instead of a tilt 
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Fig. 1: A schematic drawing of the geometry of the system. In (a) the tilt angle 9 between the normal of the 
two lamellar phases is shown. The bulk periodicity d is smaller then the local periodicity d x — dj cos(#/2) at the 
interface, y = 0. In (b) a schematic drawing of a chevron morphology with rounded V-shaped tips is shown. 

Fig. 2: Black and white contour plot of the order parameter profile, as obtained from the profile solutions of 



Sec. IV. The grayness denotes the value of the order parameter (A/B relative volume fraction). Black domains 
are A-rich, white domains are B-rich. The value of the interaction parameter is x = 1 and the average difference 
of monomer volume fractions is </>o = 0. For the small angle (9 = 20°, top plot) chevron morphology appears, 
while for 9 = 130° (bottom plot) omega takes over. 

Fig. 3: Same as Fig. 2, but with x = 0.76. Top plot is for 9 = 20° and the bottom plot is for 9 = 130°. 
The system is closer to the ODT than that of Fig. 2, and the modulations of the lamellar thickness are more 
prominent. 

Fig. 4: Black, white and gray contour plot, showing A-rich, B-rich and interfacial regions, respectively. Black 
regions are A-rich (</>q > 0.2), white regions are B-rich ((j)Q < 0.2), while gray regions marks the interfacial 
region (—0.2 < <pn < 0.2). Top plot is for 9 = 20° and in the bottom plot 9 — 130°. The interaction is set to 

x = i- 

Fig. 5: As in Fig. 4, but with \ = 0.76. Top plot is for 9 = 20° and the bottom plot is for 9 = 130°. Comparison 
with Fig. 4 shows that for fixed tilt angle 9, the omega structure and the modulations of the lamellar thickness 
are more evident closer to the ODT. 

Fig. 6: Plot of the equi-0 contour lines corresponding to <f> — 0.5, with the same conditions as in Fig. 4b. 
Although the profile shown in Fig. 4 is symmetric with respect to the interchange of A and B, clearly the 
leading and trailing edges of the lamellae are not identical. 

Fig. 7: Plot of the grain-boundary energy 7 as function of the tilt angle 9, for \ — 0.76. 9 ranges from to 



180°. The curve marked with rectangles shows the results obtained in Sec. Ill for the chevron morphology, with 
7 = f^gfZy- The curve with circles shows the results of Sec. IV. Notice that the two expressions are very similar, 



but at some intermediate values of 9 (9 ~ 50°) there is a crossover from the chevron to the omega morphology. 
Evidently, for large tilt angles the omegas cost less energy than the chevron. 

Fig. 8: In (a) the grain-boundary energy 7 of Fig. 6, is plotted as 7 1 / 3 against the tilt angle 9. The dashed 
curve with rectangles is straight for small angles, depicting perfect 7 ~ 9 3 dependence. The curve marked with 
circles shows a power exponent 2.91. In (b) the two line tension are plotted on a log-log plot. The results of 



Sec. Ill show 9 3 scaling, while the omegas show a smaller power exponent. Possible small numerical errors in 
our scheme to find the wavevector k are important for small angles (the leftmost point being at 9 = 2°), and 
this may easily change the line tension exponent from the expected 3 to 2.91. 

Fig. 9: The imaginary and real parts of the wave-vector k in expression ( |3l|) are shown in (a), for X = 1. The 
line with rectangles corresponds to the imaginary part. Only the positive imaginary part is shown. The line 
with circles corresponds to the real part of k. Dotted line shows q y , being the characteristic length scale in the 
y-direction. Notice the splitting of the real part of k in the band gaps. An enlargement of the same curve of 
the imaginary part of k, for 17° < 9 < 50° is shown in (b). The black lines indicates integer multiples of q y in 
order to show the regions of degeneracy. Three band gaps are shown, around 18°, 26° and 42°. In these gaps 
Im(fc) = imq y for some integer m. 



11 



OQ 



► 



Si 



X 



■r 



This figure "fig2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/9905378v4 



This figure "fig3.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/9905378v4 



This figure "fig4.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/9905378v4 



This figure "fig5.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/9905378v4 



0.01 





20 40 60 80 100 120 140 160 180 

e [deg] 




Fig. 8 



This figure "fig9.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/9905378v4 



